Computer algebra-perturbation solution to a nonlinear wave equation

In this paper, the higher-order asymptotic solution to the Cauchy problem of a nonlinear wave equation is found by using a computer algebra-perturbation method. The secular terms in the solution from straightforward expansions are eliminated with the straining of characteristic, coordinates and the use of the renormalization technique, and the four-term uniformly valid solution is obtained with the symbolic computation by using a computer algebra system. The comparison of the derived asymptotic solution and the numerical solution shows that they coincide with each other for smaller ε and agree quite well for larger ε (e. g., ε=0.25) Project supported by the National Natural Science Foundation of China and Shanghai Municiple Natural Science Foundation     

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

落柱式高压粘度装置的研制与压粘试验结果

作者研制了一台可测试润滑油在高压下的粘度和密度的落柱式高压粘度装置。本装置采用了新型高压粘性动密封和新型高压流变平垫静密封,经1200MPa试压、1000MPa压力下48小时保压和使用试验证明均无泄漏。用本装置所获得的压粘数据的重现性好、灵敏度高、可比性强。     

Characterization of particle-laden, confined swirling flows by phase-doppler anemometry and numerical calculation

The particle dispersion characteristics in a confined swirling flow with a swirl number of approx. 0.5 were studied in detail by performing measurements using phase-Doppler anemometry (PDA) and numerical predictions. A mixture of gas and particles was injected without swirl into the test section, while the swirling airstream was provided through a co-flowing annular inlet. Two cases with different primary jet exit velocities were considered. For these flow conditions, a closed central recirculation bubble was established just downstream of the inlet.

The PDA measurements allowed the correlation between particle size and velocity to be obtained and also the spatial change in the particle size distribution throughout the flow field. For these results, the behaviour of different size classes in the entire particle size spectrum, ranging from about 15 to 80 μm, could be studied, and the response of the particles to the mean flow and the gas turbulence could be characterized. Due to the response characteristics of particles with different diameters to the mean flow and the flow turbulence, a considerable separation of the particles was observed which resulted in a streamwise increase in the particle mean number diameter in the core region of the central recirculation bubble. For the lower particle inlet velocity (i.e. low primary jet exit velocity), this effect is more pronounced, since here the particles have more time to respond to the flow reversal and the swirl velocity component. This also gave a higher mass of recirculating particle material.

The numerical predictions of the gas flow were performed by solving the time-averaged Navier-Stokes equations in connection with the well known kε turbulence model. Although this turbulence model is based on the assumption of isotropic turbulence, the agreement of the calculated mean velocity profiles compared to the measured gas velocities is very good. The gas-phase turbulent kinetic energy, however, is considerably underpredicted in the initial mixing region. The particle dispersion characteristics were calculated by using the Lagrangian approach, where the influence of the particulate phase on the gas flow could be neglected, since only very low mass loadings were considered. The calculated results for the particle mean velocity and the mass flux are also in good agreement with the experiments. Furthermore, the change in the particle mean diameter throughout the flow field was predicted approximately, which shows that the applied simple stochastic dispersion model also gives good results for such very complex flows. The variation of the gas and particle velocity in the primary inlet had a considerable impact on the particle dispersion behaviour in the swirling flow and the particle residence time in the central recirculation bubble, which could be determined from the numerical calculations. For the lower particle inlet velocity, the maximum particle size-dependence residence time within the recirculation region was considerably shifted towards larger particles.     

复合材料光弹性分析的近似方法

本文提出了一种用于光弹性复合材料的简化应变——光学定律。按照这一简化定律。模型材料的主应变差和主应变方向只要利用光弹性实验测出的等差线与等倾线即可求得。些是一种正交异性光弹性分析的近似方法,这一方法所得结果与实验数据比较,最大误差在10％左右。由于采用简化应变——光学定律使得正交异性光弹性分析工作大为简便,因此它是一种适合于工程应用的近似方法。     

Experimental studies of suspended ferromagnetic fibers in a magnetic field

The analysis of the rotation of a ferromagnetic ellipsoid suspended in a Newtonian fluid and subjected to a uniform magnetic field is extended to include a long, slender cylindrical fiber which is magnetically saturated. Experimental observations of rotating nickel cylinders with aspect ratiosL/D ranging from 5 to 40 agree with the theoretical predictions that: (1) the proper magnetoviscous time constant for the motion is _MV = _s/µ ₀ M _s ² , (2) larger fiber aspect ratios result in considerably longer orientation times; and (3) the strength of the applied external field has only a slight effect on the overall fiber rotation, and has no effect on the maximum angular velocity achieved. Quantitative agreement of theory and experiments is obtained for fibers withL/D 20; for the shorter fibers, the theory tends to overpredict the fiber rotation rate by as much as 30%. D diameter of the cylinder - D ^P (r) position-dependent demagnetization tensor, implicitly defined in eq. (2.5) - D _xx,D _yy,D _zz volume-averaged demagnetizing factors for an ellipsoid equivalent to a uniformly magnetized cylinder, defined in eq. (2.6) - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - k geometric parameter in the hydrodynamic resistance of a body rotating in a Newtonian fluid, eq. (2.2) - L length of the cylinder - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component ofL ^(h) on the cylinder - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.4); thez-component ofL ^(m) on the cylinder - M the magnetization of a magnetic material - M _s the saturation magnitude ofM, approached by all ferromagnetic materials asH _i becomes large - r position vector of a point within a ferromagnetic body - V volume of a magnetic particle - x, y, z rectangular coordinate axes fixed in the cylinder according to figure 1 - angle of inclination of the axis of the cylinder with respect toH ₀ - shear rate - small parameter of slender body theory,=1/ln (2L/D) - _s constant viscosity of the suspending fluid - µ ₀ the magnetic permeability of free space,µ ₀=4 · 10^–7 H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - ₁ the first normal-stress coefficient - ; _z angular velocity of a rotating body; angular velocity of a cylinder about thez-axis, _z =– d/dt     

Rheological behaviour of two-phase flow of solid particles in a gas

The present investigation was concerned with the rheological behaviour of dilute suspensions of solid particles in a gas in a vertical cocurrent flow moving upwards. Starting from the experimentally determined dependence of the pressure drop on the concentration of solid particles and the Reynolds number of the carrier medium in the steady flow region, the rheological parameters were estimated using pseudo-shear diagrams. Air was the carrier medium and the dispersed phase was one of six fractions of polypropylene powder and five fractions of glass ballotini. The results show that the investigated two-phase systems have pseudoplastic character which becomes more pronounced with increases in concentration, equivalent diameter and density of solid particles in the flowing suspension. C _d coefficient of particle resistance - d _e equivalent diameter of particles - D column diameter - Fr Froude number - g gravitational acceleration - K rheological parameter - L length - n rheological parameter - p _t pressure drop due to friction - p _m total pressure drop - p _ag pressure drop due to acceleration of the gas phase - p _as pressure drop due to acceleration of the solid phase - p _g hydrostatic pressure of the gas phase - p _s specific effective weight of the dispersed phase - r radius - Re Reynolds number - Re _p Reynolds number of a particle - Re _G generalized Reynolds number - Re _G1 generalized Reynolds number relating to the end of the laminar flow region - Re _G2 generalized Reynolds number relating to the beginning of the turbulent flow region - w _z axial component of velocity - u _t steady free-fall velocity of a single particle - w average velocity - w _g average velocity of the gas phase - w _s average velocity of the dispersed phase of solid particles - relative mass fraction of solid particles - x _s volume fraction of solid particles - _g coefficient of pressure drop due to friction - µ dynamic viscosity - _g density of the gas phase - _m density of the suspension - _s density of solid particles - _ds density of the dispersed phase - _w shear stress at the wall     

基于CUDA的有限元矩阵并行装配算法研究

构建航天飞行器的结构有限元模型是准确模拟飞行仿真、完成飞行器在轨飞行阶段结构故障监测和诊断的基础。采用细长体飞行器简化梁模型，提出新的基于CUDA（Compute Unified Device Architecture）的有限元单元刚度矩阵生成和总刚度矩阵组装算法。依据梁单元矩阵的对称性，结合GPU硬件架构提出并行生成算法并进行改进。为有效减少装配时间，在装配过程中采用着色算法，提出了基于GPU（Graphics Processing Unit）共享内存的非零项组装策略，通过在不同计算平台下算例对比，验证了新算法的快速性。数值算例表明，本文算法的求解效率较高，针对一定计算规模内的模型可满足快速计算与诊断的实时性要求。     

用间接边界元方法求解不可压缩孤立翼型绕流的一种库塔条件处理方法

本文对不可压二维翼型势流绕流的边界元法求解作了分析,在对基本方程与边界条件进行数值离散化时,将库塔条件代入基本方程。按本方法编制的计算机程序对若干算例进行了验算。结果表明,本文提出的方法是可靠的,该法计算简捷、方便,占用计算机内存少,具有实用意义。     

金属材料的激光表面改性及其在摩擦学中的应用

作者对激光表面处理的各种方法和现有技术水平进行了综述介绍,并就其在改善金属材料表面摩擦磨损性能方面的研究和应用,以及今后的发展都作了简要的分析和讨论。激光相变硬化能使处理的工件形成具有表面压应力的硬质马氏体表面,因而可以降低表面的磨损速率。在金属工件上通过激光包覆高抗磨、抗热、抗腐蚀和抗疲劳的高硬质合金甚至致密陶瓷,在严格控制工件材料对包覆合金稀释的情况下,可使底材获得高的摩擦学性能。激光合金化通过选择合金元素和基底材料能有选择地改进低成本工件的表面,使其具有优异的物理、化学和机械性能。激光上釉是通过细化铸造组织、减少偏析及形成高度过饱和固溶体等亚稳定相乃至非晶态而提高了材料的耐磨性、抗氧化性和抗腐蚀性能。激光冲击硬化是通过在材料表层内形成冲击波而引起“损伤”来改变表层的组织和性能。